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Modeling  of  Plasma-Induced  Ignition  and 
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Iain  D.  Boyd 

Department  of  Aerospace  Engineering.  University  of  Michigan 


Abstract —  Much  progress  has  been  made  in  the  work  performed 
in  this  grant  toward  the  development  of  an  end-to-end  computer 
model  of  an  electrothermal  chemical  gun.  Phenomena  that  must 
be  considered  in  an  electrothermal  chemical  gun  model  include 
the  initial  capillary  plasma  properties,  the  plasma-air  interaction, 
plasma  sheath  effects,  and  the  plasma-propellant  interaction  itself. 
This  report  presents  an  overview  of  recent  progress  made  modeling 
these  various  phenomena.  Also  discussed  are  methods  under 
consideration  to  incorporate  surface  chemistry  effects  into  the 
plasma-propellant  interaction  model. 

Index  Terms — capillary,  ETC,  heat  flux,  modeling,  plasma, 
plasma  chemistry,  plasma  sheath. 


I.  Introduction 

There  has  been  much  interest  in  recent  years  in  the  design 
and  implementation  of  electrothermal  chemical  (ETC)  guns. 
An  ETC  gun  is  a  solid  propellant  based  artillery  piece  in 
which  the  conventional  ignition  system  is  replaced  with  a 
capillary  plasma  source.  Among  the  enhancements  encountered 
in  plasma  based  ignition  systems  are  reduced  ignition  delay 
time,  highly  repeatable  ignition  time,  and  enhanced  burning 
and  combustion  of  the  solid  propellant  [1],  [2], 

Work  is  continuing  in  the  development  of  an  end-to-end  com¬ 
puter  model  of  an  ETC  gun,  focused  primarily  on  modeling  the 
physics  of  the  plasma-propellant  interaction.  To  date,  models 
have  been  developed  that  address  the  capillary  plasma  source 
[3],  [4],  the  plasma-air  chemistry  of  the  expanding  capillary 
plasma  jet  into  the  combustion  chamber  [5],  the  plasma- 
propellant  interaction  via  a  coupled  ablation  and  thermal  model 
[6],  and  the  convective  heat  flux  to  the  propellant  bed  by  means 
of  a  collisional  plasma  sheath  model  [7],  This  effort  represents 
the  first  attempt  to  create  a  model  of  the  ETC  phenomena, 
from  the  initial  capillary  firing  through  the  plasma-propellant 
interaction  leading  up  to  propellant  ignition. 

In  this  report,  we  present  an  overview  of  the  work  done  on 
the  ETC  model,  including  the  major  conclusions  of  the  efforts. 
Section  II  addresses  the  capillary  model.  Section  III  summarizes 
CFD  studies  of  the  plasma-air  interaction.  The  previously 
published  plasma-propellant  interaction  model  is  covered  in 


Section  IV,  which  is  loosely  coupled  to  a  collisional  plasma 
sheath  model  described  in  Section  V.  A  summary  of  current 
work  on  propellant  surface  chemistry  is  included  in  Section 
VI,  with  discussion  on  final  model  integration  in  Section  VII. 
Concluding  remarks  are  presented  in  Section  VIII. 


II.  Capillary  Model 

The  capillary  model  [3],  [4]  uses  kinetic  ablation  theory 
for  the  wall  material  [8],  [9],  coupled  to  a  non-equilibrium 
hydrodynamic  layer.  Solving  the  energy,  momentum  and  mass 
conservation  equations  in  the  hydrodynamic  layer  allows  for 
the  determination  of  the  plasma  density  and  temperature  at 
the  polyethylene  surface.  The  kinetic  ablation  model  then 
determines  the  ablation  rate,  which  in  turn  results  in  updated 
plasma  density,  pressure,  and  electron  and  surface  temperatures. 

The  model  is  for  a  one  dimensional  axisymmetic  capillary 
of  length  L  and  internal  radius  R.  The  spatial  variable  x  varies 
from  the  back  of  the  capillary  to  the  exit  plane,  at  x  =  L. 
The  model  assumes  that  all  flow  parameters  are  uniform  along 
the  radial  direction  in  each  layer.  The  energy  equation  in  the 
hydrodynamic  layer  is 

/  Os  Os  \  0U 

p{oi  +  u^)  =  -p^  +  Q^Q^QF-  (1) 

Fluid  density,  velocity,  and  pressure  are  given  by  p,  U  and  P, 
respectively.  In  (1),  e  =  (3/2 ){T/m)  +  (U2/ 2),  where  T  is 
the  plasma  temperature,  and  to  is  the  mass  of  an  average  fluid 
particle.  The  influx  of  energy  to  the  plasma  due  to  Joule  heating 
is  represented  by  Qj.  Energy  losses  due  to  plasma  radiation  and 
convection  to  the  wall  are  given  by  Qr  and  Q  /--,  respectively. 
It  is  assumed,  due  to  the  high  pressure  and  small  scale  of  the 
capillary,  that  pressure  is  constant  in  the  hydrodynamic  layer, 
thus  OP /Ox  =  0.  Previous  simulations  and  theoretical  estimates 
indicate  that  the  plasma  temperature  varies  only  slightly  in  the 
capillary,  thus  it  is  assumed  that  OT/Ox  =  0.  [3] 

Solving  the  energy,  momentum  and  mass  conservation  equa¬ 
tions  in  the  hydrodynamic  layer  yields  the  following  relation 


for  the  ablation  rate,  F, 
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where  k  is  the  Boltzmann  constant  and  n  is  the  fluid  number 
density.  Subscripts  1  and  2  represent  values  at  the  Knud- 
sen/hydrodynamic  layer  and  the  hydrodynamic/quasi-neutral 
plasma  layer  boundaries,  respectively  [3].  Note  that  (2)  is  only 
real  if  nikTi  >  U2kT2  and  n\  >  712.  If  these  conditions 
are  not  met  the  implication  is  that  the  backflux  is  higher  than 
the  flux  from  the  wall  into  the  bulk  plasma.  When  this  is  the 
case,  deposition  to  the  polyethylene  surface  takes  place,  and  the 
total  ablated  mass  decreases.  Both  conditions  are  met  during 
the  discharge,  thus  the  capillary  model  accommodates  both 
situations.  Polyethylene,  and  therefore  the  capillary  plasma,  is 
composed  of  carbon  and  hydrogen.  Studies  of  film  deposition 
have  shown  that  carbon  is  likely  to  deposit  and  accumulate  on 
the  surface,  while  hydrogen  is  not.  To  account  for  this  the  model 
includes  a  parameter  v ,  which  is  the  fraction  of  the  backflux  that 
deposits  on  the  surface,  with  v  =  1  being  the  most  deposition, 
and  v  =  0  neglecting  deposition  entirely.  Previous  simulations 
and  experiments  found  good  agreement  with  v  =  0.6  [4]. 


As  inputs  the  model  requires  capillary  length  and  internal 
radius,  a  value  for  the  backflux  parameter  v ,  and  a  time  varying 
electric  current  profile.  Polynomial  curve  fits  to  experimentally 
measured  discharge  currents  can  be  input  into  the  model  to 
attempt  to  match  specific  experimental  data. 


As  outputs  the  capillary  model  generates  time  varying  data 
for  the  capillary  outflow  density,  pressure,  temperature,  species 
composition,  and  total  capillary  ablated  mass.  Outflow  velocity 
is  assumed  to  be  sonic  as  a  boundary  condition  in  the  hy¬ 
drodynamic  model.  The  profile  of  the  outflow  pressure  closely 
follows  that  of  the  discharge  current.  Fig.  1  shows  a  polynomial 
fit  to  an  experimental  current  profile  and  the  resulting  outflow 
pressure  and  plasma  temperature,  Te,  over  a  220  /is  discharge. 
The  primary  accomplishment  of  this  model  is  it’s  ability  to 
accurately  reproduce  or  predict  the  total  ablated  mass,  and  time 
varying  outflow  pressure  and  species  composition. 


III.  Plasma-Air  Chemistry 
During  the  capillary  discharge,  the  generated  plasma  exits  the 
capillary  nozzle  and  either  comes  into  immediate  contact  with 
the  propellant,  or  expands  in  an  open  air  filled  chamber  before 
encountering  the  propellant  bed,  depending  on  the  propellant 
packing  option  being  used  [2],  [10],  Our  simulations  indicate 
that,  for  typical  ETC  conditions,  if  the  plasma  jet  is  allowed  to 
expand  into  air,  plasma-air  chemistry  will  occur  if  the  capillary- 
propellant  distance  is  less  than  approximately  10  mm  [5].  For 
geometries  with  plasma-air  interaction  regions  on  this  scale, 
plasma-air  chemistry  should  be  taken  into  account. 

A  plasma-air  chemistry  model  based  on  previous  experimen¬ 
tal  and  theoretical  work  was  developed  that  includes  20  species 


Fig.  1.  Example  input  and  output  (dashed  lines)  from  the  capillary  model, 
taken  from  [5]. 
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Fig.  2.  Mole  fraction  of  OH  in  the  plasma  jet,  taken  from  [5].  Production  is 
limited  to  a  thin  plasma-air  reaction  region  at  the  periphery  of  the  expanding 
jet. 


and  41  chemical  reactions  [5].  This  model  was  used  in  a  chem¬ 
ically  reacting  CFD  code  developed  at  the  University  of  Michi¬ 
gan,  known  as  LeMANS.  LeMANS  is  a  parallelized,  implicit 
Navier-Stokes  CFD  code  that  uses  a  modified  Steger- Warming 
Flux  Vector  Splitting  scheme  for  calculation  of  inviscid  fluxes 
between  mesh  volumes,  is  second  order  accurate  in  space,  with 
time  integration  performed  using  a  point  implicit  method  [11], 
[12],  The  species  included  in  the  chemistry  model  are  N2,  O2, 
NO,  N,  O,  OH,  CO,  CO2,  H20,  CH,  HCO,  NH,  H2,  C,  H, 
H+,  C+,  N+,  0+,  and  electrons.  Electron  impact  ionization  is 
assumed.  The  output  from  the  capillary  model  (Section  II)  is 
used  as  time  varying  input  to  LeMANS.  Electron  mobility  is 
limited  in  the  model  to  ensure  quasi-neutrality  is  maintained 
everywhere  in  the  simulation.  Fig.  2  shows  a  representative 
plasma-air  chemistry  result  from  the  axisymmetric  plasma  jet 
simulation.  We  can  see  that  the  production  of  OH  occurs  in 
a  thin  region  where  the  hydrogen  of  the  expanding  plasma 
jet  encounters  the  oxygen  in  the  ambient  air,  as  predicted  by 
theory.  As  the  plasma  jet  expands,  this  reaction  region  is  pushed 
outward.  The  sickle-like  shape  near  X=20  mm,  Y=10  mm  is  a 
recirculation  region  [5]. 


JA2  no  (m 3) 


As  part  of  the  plasma-air  study,  a  bayonet  tube  was  simulated 
that  approximated  the  physical  dimensions  used  in  [10],  A 
bayonet  tube  is  inserted  into  a  cylindrical  bore  in  a  propellant 
block.  The  tube  has  holes  spaced  along  it’s  length  which 
provide  controlled  ignition  sites  and  allows  ignition  to  occur 
not  only  at  the  end  of  the  propellant  closest  to  the  capillary, 
but  at  the  far  end  as  well.  Data  from  the  simulated  “holes” 
clearly  shows  that  by  the  time  the  plasma  wave  reaches  the 
end  of  the  tube  a  significant  amount  of  plasma-air  chemistry 
has  occurred.  This  indicates  that  in  a  bayonet  tube-like  packing 
configuration  the  ignition  mechanism  at  one  end  of  the  tube 
may  be  different  from  the  mechanism  at  the  other  end,  due  to 
the  drastically  different  plasma  chemical  composition  [5], 
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One  of  the  main  benefits  of  the  plasma-air  CFD  studies  is  the 
ability  to  predict  the  plasma  jet  chemical  composition  under  dif¬ 
ferent  geometries  and  capillary  inflow  conditions.  Experimental 
attempts  to  measure  the  composition  of  the  plasma  jet  have 
yielded  few  promising  results  [5],  The  results  of  the  simulations 
clearly  indicate  that  small  capillary-propellant  distances  allow 
us  to  entirely  neglect  plasma-air  chemistry,  which  yields  insight 
when  choosing  propellant  packing  options  for  practical  ETC 
application. 


IV.  Plasma-Propellant  Interaction  Model 
The  primary  focus  of  the  ETC  model  is  the  simulation  of 
the  plasma-propellant  interaction  (PPI).  This  is  accomplished 
via  a  kinetic  theory  based  ablation  model  adapted  from  work 
by  Keidar  et  al  [9],  To  model  a  propellant  we  require  the 
vapor  pressure,  enthalpy  of  sublimation  (AH),  specific  heat 
at  constant  pressure  ( Cp ),  the  thermal  conductivity  (A)  and 
diffusivity  (a).  Simulations  [6]  were  run  comparing  two  propel¬ 
lants  under  consideration  for  ETC  application.  JA2  is  a  double¬ 
base  propellant,  which  is  optically  semi-transparent,  allowing 
radiation  to  penetrate  in-depth.  XM39  is  a  nitramine  composite 
propellant,  is  optically  opaque  to  most  wavelengths,  and  does 
not  allow  significant  radiation  penetration. 

For  the  PPI  study,  a  range  of  plausible  bulk  plasma  densities 
are  chosen.  Plasma  number  density  ( n0 )  ranged  from  1021-1024 
m-3.  The  ablation  model  is  similar  to  that  used  in  the  capillary 
model  (Section  II),  with  the  ablation  rate  being  a  function 
of  bulk  plasma  density  and  propellant  surface  temperature. 
Surface  temperature  is  obtained  via  a  one-dimensional  thermal 
model  representing  the  propellant  bed.  Experimental  data  was 
available  for  total  ablated  mass  of  a  small  propellant  disc  after 
exposure  to  a  capillary  plasma  jet,  this  experimental  geometry 
was  used  in  the  thermal  model  [13].  The  propellant  disc  was 
of  thickness  L  =  4  mm  and  was  positioned  in  a  sample  holder 
taken  to  be  a  thermal  sink.  It  is  assumed  that  the  experiment  was 
conducted  at  one  atmosphere  pressure  and  initial  temperature 
of  T0  =  298  K.  The  thermal  model  and  boundary  conditions 
are 

dtT(x,t)  =  adxT(x,t), 

T(0,t)  =  T0,  T(x,0)  =  T0,  (3) 


Fig.  3.  Plasma-propellant  interaction  results  for  ablated  mass  versus  propellant 
surface  temperature  for  JA2,  taken  from  [6]. 


4.25  r 

4  7 

"e 

3.75  - 

CO 

- 

b 

T— 

3.5  - 

O- 

X 

3 

3.25  - 

LL 

3  7 

+-• 

CO 

0) 

X 

2.75  - 

2.5 


10" 


JA2 

XM39 


10"' 


10"' 


Plasma  Density,  n0  (m‘ ) 


10" 


Fig.  4.  Constant  surface  heat  flux  necessary  to  reproduce  the  experimental 
ablated  mass  as  a  function  of  bulk  plasma  density  [6], 


dxT{L,t)  =  -X -\q  -  AHT  -  CP(TS  -  T0)T), 

where  Ts  =  T(L,t),  and  q  is  the  surface  heat  flux  to  the 
propellant,  which  is  assumed  constant  in  time  and  space.  The 
objective  of  the  model  is  the  determination  of  the  effective 
surface  heat  flux  over  the  course  of  the  discharge,  q.  The  re¬ 
sulting  data  can  be  presented  as  plots  of  total  propellant  ablated 
mass  as  a  function  of  peak  propellant  surface  temperature,  at 
different  bulk  plasma  densities,  as  in  Fig.  3.  Note  that  there  are 
multiple  solutions  for  each  plasma  density  in  terms  of  the  peak 
surface  temperature.  Knowing  the  experimentally  determined 
ablated  mass  allows  us  to  determine  the  necessary  heat  flux  as 
a  function  of  plasma  density.  The  experimental  ablated  mass  for 
JA2  of  5.3  mg  is  presented  as  a  horizontal  line  in  Fig.  3.  Fig. 
4  presents  data  from  the  model  that  matches  the  experimental 
ablated  masses  for  JA2  and  XM39,  representing  the  constant 
surface  heat  flux  necessary  to  match  the  experimental  results. 
We  consider  the  ability  to  determine  the  exact  bulk  plasma 
density  when  we  discuss  model  integration  in  Section  VII. 

The  plasma-propellant  interaction  model  led  to  several  in¬ 
teresting  conclusions.  Fig.  4  clearly  shows  that  JA2  will  con- 


sistently  have  a  higher  heat  flux  than  XM39.  Fig.  3  and  the 
corresponding  plot  for  XM39  show  that  XM39  has  a  higher 
surface  temperature  than  JA2  at  a  fixed  plasma  density  [6]. 
The  downward  trend  in  Fig.  4  is  caused  by  the  fact  that  as 
plasma  density  increases,  the  heat  flux  must  decrease  in  order 
to  match  the  experimental  ablated  mass,  clearly  illustrated  by 
examining  ablation  rate  contours  [6].  Of  significant  interest 
is  the  difference  between  the  two  heat  fluxes  in  Fig.  4.  If 
the  capillary  generates  an  identical  plasma  at  each  firing  (a 
valid  assumption)  then  why  is  the  total  heat  flux  different? 
The  most  plausible  explanation  involves  the  propellant’s  optical 
properties.  JA2’s  ability  to  access  plasma  radiation  as  a  heat 
source  allows  it  to  tap  into  more  of  the  available  thermal 
energy  of  the  plasma,  while  XM39’s  opacity  reflects  most  of 
this  radiative  energy  back  into  the  plasma.  It  is  also  possible 
that  plasma  surface  chemistry  is  responsible  for  part  of  this 
difference,  and  is  being  explored  currently,  as  we  will  discuss 
in  Section  VI. 


V.  Collisional  Plasma  Sheath  Model 

The  work  modeling  the  plasma-propellant  interaction  yielded 
the  total  heat  flux  incident  on  the  propellant  surface  [6],  In  order 
to  determine  what  phenomena  contribute  to  the  total  heat  flux, 
a  collisional  plasma  sheath  model  was  developed  to  find  the 
convective  heat  flux  to  the  propellant  bed  [7], 

A  plasma  sheath  develops  at  any  plasma-solid  interface.  High 
speed  electrons  strike  the  propellant  before  the  slower,  heavier 
ions,  causing  a  negative  potential  to  build  up  on  the  propellant. 
This  negative  potential  grows  until  the  electron  and  ion  fluxes 
balance,  and  a  floating  potential  is  reached.  This  plasma  sheath 
causes  ions  to  accelerate  as  they  approach  the  wall  due  to 
the  potential  field,  thus  increasing  their  kinetic  energy.  Ion- 
neutral  collisions  in  the  sheath  decrease  the  ion  drift  velocity, 
decreasing  the  ion  kinetic  energy.  These  two  factors  determine 
the  convective  ion  flux  to  the  propellant  bed  in  ETC  application. 

The  one-dimensional  sheath  model  assumes  Boltzmann  elec¬ 
trons,  a  single  fluid  ion  species,  and  a  uniform  neutral  back¬ 
ground  density.  The  neutral  density  is  found  from  the  PPI 
model,  where  the  neutral  number  density  is  given  by  nn  = 
Pv(Ts)/(kTs),  where  Py(Ts)  is  the  propellant  vapor  pressure 
at  surface  temperature  Ts,  as  determined  by  the  PPI  thermal 
model  [6],  In  the  results  presented  here,  a  bulk  plasma  temper¬ 
ature  of  Te  =  1.5  eV  is  assumed.  The  only  remaining  unknown 
is  the  velocity  of  the  ions  as  they  enter  the  sheath.  Theoretical 
limits  for  a  collisional  plasma  sheath  indicate  that  the  inflow 
velocity  should  be  less  than  or  equal  to  the  Bohm  velocity, 
given  as  Cs  =  \J eTe/nii,  where  rnt  is  the  ion  mass.  We  non- 
dimensionalize  the  inflow  velocity  as  Ui0  =  Vi0/Cs,  and  select 
representative  values  of  Ui0  =0.3,  0.6,  and  1.0.  The  convective 
heat  flux  to  the  propellant  bed  is  given  as 

Qconv  (tr0C s^io)  ^  2  ( 6?s ‘Uj:W )  “b  2c7).  eTer/wJ  ,  (4) 

where  n0  is  the  bulk  plasma  number  density,  and  ?/„;  is  the 
non-dimensional  floating  potential,  calculated  via  an  analytical 


Fig.  5.  Convective  heat  flux  to  XM39  for  n0  =  3  X  1022m  3  for  different 
ion  inflow  velocities,  taken  from  [7], 


expression  [7].  The  only  unknown  in  (4)  is  UiW,  the  non- 
dimensional  ion-wall  impact  velocity,  which  we  obtain  from 
the  sheath  solution. 

Fig.  5  shows  the  convective  heat  flux  to  XM39  for  a  bulk 
plasma  density  of  n0  =  3  x  1022?n~3  and  different  ion  inflow 
velocities,  where  At  =  1  is  the  length  of  the  plasma  pulse  from 
the  PPI  work  of  Section  IV.  We  can  see  that  the  convective 
heat  flux  changes  curvature  in  Fig.  5,  and  appears  to  approach 
a  limit  near  At  =  0,  and  clearly  reaching  a  different  limit  well 
before  At  =  1.  The  two  extremes  are  the  collisionless  and  fully 
collisional  plasma  sheath  limits.  Analytical  expressions  exist 
for  the  ion-wall  impact  velocities  for  both  fully  collisional  and 
collisionless  plasma  sheaths;  these  analytical  limits  are  shown 
in  Fig.  6  for  the  uto  =  1  case.  During  the  early  portion  of 
the  discharge,  the  propellant  surface  temperature,  and  thus  the 
neutral  density,  is  low,  resulting  in  a  low  collisionality  in  the 
sheath.  By  the  end  of  the  pulse  the  neutral  density  is  high 
enough  that  the  sheath  is  fully  collisional.  However,  the  sheath 
remains  in  a  transition  regime  for  most  of  the  plasma  pulse. 
Since  JA2  has  a  higher  ablated  mass,  results  for  that  propellant 
typically  do  not  come  as  close  to  the  collisionless  limit  as  do 
the  results  for  XM39. 

The  main  conclusions  from  the  collisional  plasma  sheath 
model  include  the  fact  that  XM39  receives  a  higher  convective 
heat  flux  than  JA2  for  a  given  plasma  density  and  ion  inflow 
velocity.  Given  that  the  PPI  model  shows  that  JA2  has  a  higher 
total  heat  flux,  this  means  that  JA2  receives  enough  plasma 
radiation  to  overcome  XM39’s  higher  convective  heat  flux, 
indicating  that  plasma  radiation  is  a  significant  heat  transport 
mechanism  in  the  PPI.  As  in  the  PPI  model,  the  bulk  plasma 
density,  n0,  is  treated  as  a  parameter  as  it  is  not  known  at  this 
time.  We  find  that,  depending  on  the  value  of  the  bulk  plasma 
density,  convective  heat  flux  can  represent  a  significant  portion 
of  the  total  heat  flux  to  the  propellant  bed  [7]. 
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Fig.  6.  Comparison  of  simulation  result  to  limits  for  collisionless  and  fully 
collisional  ion- wall  impact  velocity,  taken  from  [7], 


Fig.  7.  Vapor  pressure  of  JA2.  XM39,  and  carbon  as  a  function  of  surface 
temperature  [6],  [15].  Note  that  the  carbon  vapor  pressure  is  negligible  when 
compared  to  that  of  the  propellant. 


VI.  Surface  Chemistry 

Current  work  on  the  ETC  model  is  focused  on  adding  surface 
chemistry  effects  to  the  PPI  model  and  final  integration  of  all 
the  submodels  (Sections  II- V).  The  focus  of  the  investigation 
into  surface  chemistry  phenomena  is  carbon  film  deposition  on 
the  propellant  surface.  The  propellants  are  composed  primarily 
of  carbon,  hydrogen,  oxygen  and  nitrogen.  Experimental  abla¬ 
tion  studies  frequently  show  carbon  deposition  to  the  surface  of 
the  ablated  material  [14],  Carbon  deposition  to  the  propellant 
surface  will  affect  the  ablation  rate  in  two  ways.  First,  the 
ablation  of  propellant  will  decrease,  as  the  carbon  film  forms 
a  patchy  covering  on  the  propellant  which  shields  exposure  to 
the  plasma.  Secondly,  as  carbon  attaches  itself  to  the  propellant, 
or  to  the  already  existing  carbon  film,  a  new  chemical  bond 
is  formed  which  extracts  energy  from  the  material  surface. 
This  energy  loss  can  be  accounted  for  in  the  thermal  boundary 
condition  in  (3).  The  energy  lost  in  this  way  will  decrease  the 
propellant  surface  temperature,  causing  a  further  decrease  in 
the  ablation  rate. 

As  a  first  step  toward  implementing  film  deposition  in  the  PPI 
model,  the  ablation  rate  for  carbon  was  found  via  the  ablation 
model  used  in  Sections  II  and  IV,  which  is  based  on  the  material 
vapor  pressure.  It  was  found  for  the  surface  temperature  range 
predicted  by  the  PPI  model  (300  K  <  Ts  <  1200  K)  and  for 
the  conditions  of  the  experiment  used  to  fix  the  total  ablated 
mass  [6],  [13],  that  carbon  ablation  is  negligible.  This  is  due 
to  carbon’s  low  vapor  pressure  in  the  temperature  range.  Fig. 
7  shows  the  vapor  pressure  of  JA2  and  XM39,  taken  from 
[6],  compared  to  that  of  carbon,  taken  from  [15],  over  the 
temperature  range  of  interest.  We  see  that  the  carbon  vapor 
pressure  is  much  less  than  the  vapor  pressure  of  either  JA2  or 
XM39,  resulting  in  negligible  carbon  ablation. 

The  rate  of  carbon  film  growth  is  thus  the  critical  aspect 
of  the  model,  as  any  carbon  that  is  deposited  will  not  re-ablate 
during  the  PPI.  Typically,  the  amount  of  material  deposited  back 
onto  the  ablating  surface  is  calculated  using  a  fitting  parameter 
[4],  The  ability  to  use  a  fitting  parameter  does  not  present  itself 


here,  thus  we  must  model  this  in  some  way.  One  approach 
being  considered  is  to  use  the  carbon-carbon  bonding  energies 
and  rates  to  estimate  a  probability  that  a  carbon  atom  will 
attach  to  the  propellant  surface,  forming  a  film.  These  efforts 
are  currently  underway. 


VII.  Model  Integration 

Here  we  discuss  the  integration  of  all  submodels  into  an 
end-to-end  ETC  model.  As  we  have  seen  from  Section  IV, 
the  PPI  model  requires  the  plasma  density  at  the  propellant 
surface.  Previously,  this  information  was  not  available,  and  a 
parametric  approach  was  used,  as  shown  in  Fig.  4.  It  thus 
became  a  high  priority  to  determine  the  exact  plasma  density  at 
the  propellant  surface.  This  is  to  be  done  by  one  of  two  different 
methods,  depending  on  the  propellant  packing  option  used. 
Either  the  capillary  nozzle  is  in  direct  (or  near  direct)  contact 
with  the  propellant  surface,  in  which  case  the  capillary  output 
can  be  used  directly  to  find  the  plasma  density,  or  the  capillary- 
propellant  distance  is  greater  than  approximately  10  mm,  in 
which  case  the  plasma-air  chemistry  model  should  be  used 
to  determine  the  plasma  density,  especially  if  there  is  a  large 
chamber  which  allows  the  plasma  jet  to  undergo  expansion. 
Although  the  former  option  is  most  likely  in  practical  ETC 
application,  we  must  still  use  the  experimental  ablated  mass  to 
fix  the  total  heat  flux  in  the  PPI  model,  and  that  data  comes 
from  expanding  plasma  jet  studies  [13],  Thus  we  will  use 
the  CFD  plasma-air  chemistry  model  to  determine  the  plasma 
density  at  the  propellant  surface,  keeping  in  mind  that  using  the 
capillary  model  to  find  the  plasma  density  directly  is  simpler 
in  comparison. 

Note  that  the  PPI  results  were  utilized  in  the  convective  heat 
flux  model  (Section  V),  meaning  that  if  the  plasma  density  and 
total  heat  flux  are  known,  then  the  convective  heat  flux  will  also 
be  known.  The  surface  chemistry  model  described  in  Section 
VI  will,  through  the  energy  loss  term  in  the  thermal  model, 
decrease  the  total  heat  flux  to  the  propellant  bed,  although 


by  how  much  is  not  clear  at  this  time.  The  PPI  model  will 
readily  accommodate  a  time  varying  plasma  density,  and  using 
the  experimental  ablated  mass  we  will  obtain  a  more  accurate 
result  for  the  total  heat  flux.  The  surface  temperature  versus 
time  profile  generated  by  the  PPI  thermal  model  can  then  be 
used  in  the  convective  heat  flux  model  to  determine  exactly 
how  important  convective  heating  is  in  the  PPI  by  allowing  us 
to  calculate  the  exact  ratio  of  convective  to  total  heat  flux. 


[15]  F.  R.  Leider,  O.  H.  Krikorian,  and  D.  A.  Young,  “Thermodynamic 
Properties  of  Carbon  Up  To  The  Critical  Point,”  Carbon ,  vol.  11,  pp. 
555-563,  1973. 


VIII.  Conclusion 

Presented  was  the  progress  toward  an  end-to-end  electrother¬ 
mal  chemical  gun  model.  The  model  covers  the  generation 
of  the  plasma  in  the  capillary  plasma  source,  the  plasma- 
air  interaction  of  the  expanding  plasma  jet,  and  the  plasma- 
propellant  interaction  leading  up  to  propellant  ignition.  The 
model  yields  good  results  when  compared  to  existing  theoretical 
and  experimental  limits.  The  results  of  the  various  simulations 
yield  insight  into  the  design  of  practical  electrothermal  chemical 
guns.  Near  term  additions  to  the  model  and  integration  of  all 
submodels  are  also  discussed. 
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